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f^ Abstract 

We present a simple micromanipulation technique to transfer suspended graphene flakes onto 
any substrate and to assemble them with small localized gates into mechanical resonators. The 
mechanical motion of the graphene is detected using an electrical, radio-frequency (RF) reflection 
readout scheme -where the time- varying graphene capacitor reflects a RF carrier at / = 5 — 6 GHz 
producing modulation sidebands at /zt/m- A mechanical resonance frequency up to /„ = 178 MHz 
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H is demonstrated. We find both hardening/softening Duffing effects on different samples, and obtain 

'^ a critical amplitude of ~ 40 pm for the onset of nonlinearity in graphene mechanical resonators. 

^ Measurements of the quality factor of the mechanical resonance as a function of DC bias voltage 

Vdc indicate that dissipation due to motion-induced displacement currents in graphene electrode 

^ is important at high frequencies and large Vdc- 
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Graphene is a perfect two dimensional crystal with high Young's modulus i? ~ 1 
TPa PP and exteremly low mass, which makes it ideal for high-frequency, high-(5 nano- 
electromechanical systems (NEMS). It is a promising material for extremely low mass 
sensors ^0n4^ and for mechanical resonators at the quantum limit [5H7]. An experimental 
demonstration of a monolayer graphene resonator using optical methods [8] was achieved 
soon after the discovery of this material [9]. Compared with optical methods [HlllO], electri- 
cal detection schemes pT| IT2] are more compatible with microelectronic applications and, 
moreover, they facilitate easier studies of fundamental phenomena at low temperatures, 
where higher Q values and better sensitivity are obtainable. 

Several methods for measuring and fabricating high frequency NEMS have been developed 
O [T3142T] . They all share the challenges brought by downsizing towards the sub micron-scale, 
since the vibration amplitude under constant drive force of the mechanical bar-type resonator 
diminishes proportional to length squared [21 [22] ■ The optimization of a NEMS-resonator 
geometry is dependent essentially on the read-out scheme to be employed. When the use 
of capacitive techniques with LC matching circuits [151 116] is envisaged, in most cases, de- 
creasing the parasitic capacitance is crucial for increasing sensitivity. Unfortunately, using 
common preparation methods for suspended graphene samples, such as undercut etching 
of a sacrificial layer {e.g. Si02) [23] or random exfoliation over predefined trenches [21], 
it is hard to make small-capacitance, localized gates for graphene mechanical resonators. 
Generally speaking, the undercut etching leads to high parasitic capacitance [25], while ran- 
dom exfoliation suffers from low success rate. There have been some pioneering attempts 
on transferring carbon nanotubes [26] and graphene pieces [27H32] with polymer films, but 
realizing suspended structures after transfer remains difficult. In this work, we have devel- 
oped the polymer transfer method into a micron-scale, e-beam patterned stamp technique, 
which allows us to move individual suspended graphene flakes and assemble them with small 
localized gates into electrically-controlled mechanical resonators. We demonstrate, for the 
first time, the dispersive readout scheme for graphene mechanical resonators. 

The basic steps of our suspended stamp technique are shown in Fig. flta). First, high 
quality tape-exfoliated monolayer graphene pieces are located on a silicon substrate covered 
with 275 nm thick Si02 (confirmed with Raman spectroscopy). 50 nm gold electrodes 
are deposited on the graphene pieces by a normal e-beam lithography (EBL) and lift-off 
procedure. Then PMMA is spun on the chip again and a 2nd EBL step is applied to define 
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FIG. 1. Fabrication of suspended graphene mechanical resonators, (a) Schematic drawing of our 
suspended stamp transfer technique, (b) SEM image of a suspended graphene in a PMMA stamp; 
the scalebar equals 1 fiui. (c) Optical image of a graphene mechanical resonator on a sapphire 
substrate with localized gate electrode; the scalebar corresponds to 5 fim. 



patterns such as stamps and holes in it. Next, we peel off the whole PMMA membrane 
(~ 4mm X 4mm) from the initial substrate by etching away the Si02 in 1% HF solution, 
rinse it with DI water, and fish the whole membrane out of water using a ring-shaped frame. 
The graphene pieces and gold electrodes are both embedded in the PMMA stamps, which 
are linked with the whole PMMA membrane by some weak joints on the edges, as shown in 
Fig. [l]^a). In fact, the graphene pieces are clamped by the gold electrodes and suspended 
across the holes in the PMMA stamps, as shown in Fig. [lib). We can locate any individual 
stamp under an optical microscope. Then using a fine-tipped glass needle controlled by a 
micromanipulator [53], we can break the joints to pick up a specific stamp, transfer and 
lay it onto any desired target spot with a precision of ~ 1 /im. Although the suspended 
monolayer graphene pieces are invisible under the optical microscope, the position can be 
well judged by the visible electrodes and holes in the stamp. During the transfer, we invert 
the stamp thus bringing the graphene and electrodes from the bottom to the top of the 



stamp, so that they are supported by the PMMA and will be kept suspended over the gate 
after assembly. The typical gap distance between the gate and the suspended graphene is 
0.1-0.5 /im, which is controlled by the thickness of the PMMA. Finally, to make electrical 
contact between the hanging gold electrodes and the premade electrodes of the target circuit, 
we press the hanging electrodes through windows in the PMMA down to punch contact pads 
using the same glass needle. The punch method results in good ohmic gold-gold contacts 
in our experiments. We make the on-chip target electrodes via normal EBL using Ti/Au 
deposition on a sapphire substrate. The suspended graphene is grounded via punch contacts 
and the localized gate below it is bound to the signal line. A typical optical image of our 
sample is shown in Fig. U[c). Since the graphene pieces are suspended during the entire 
store and transfer process, the success rate depends on the size of the suspended part. At 
present, our success rate for samples below 2 fim. is about 50%. 

Cavity-based capacitive readout methods have proven very efficient in studying microme- 
chanical resonators, but their application becomes increasingly difficult when the operating 
frequency is increased towards 1 GHz [15]. Compared with mixing techniques, the loss in 
sensitivity due to larger measurement frequency is compensated by the possibility of working 
without coherent excitation for driving the mechanical oscillations. The main problem in 
high frequency readout is stray capacitance that easily masks the variation in capacitance 
SC ^ 1 fF induced by mechanical vibrations. Parasitic capacitance can be eliminated to a 
high degree by tuning the circuitry by an on-chip inductance (TB] , but even then the sensitiv- 
ity may remain moderate because of the stray capacitance of the inductor itself. A detailed 
analysis of the ideas involved in the RF capacitive readout method can be found in a recent 
review \\A\ . In this work, we have employed a vr-matching network where the Al bond wire 
forms the inductive element, in a manner shown in the inset of Fig. |2j The vr-matching 
circuit enhances the sensitivity of the measurement, but the actual result is very sensitive 
to the bond wire length and the optimization must be done by trial and error. 

For a clamped-clamped graphene mechanical resonator with length L, the time vary- 
ing average gap between the gate and the graphene can be written as D(t) = D^q + 
cos{(jJmt)\ ^Q A{y)dy when vibrating at a mechanical resonance frequency fm = Um/^n. 
Here Deq is the average gap at the equilibrium position, and A{y) denotes the amplitude of 
the flexural mode. The oscillation of the gap width gives rise to a time-varying capacitance 
Cg{t) = Ceq + 5C cos{uJmt) whcrc 5C = ^\ Jq A{y)dy- According to lumped-element mod- 
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FIG. 2. Mechanical resonance frequency fm for sample #1 measured at 4.2K. The amplitude of 
the sideband V± is given by the color scale in Volts as a function of driving frequency / and DC 
gate voltage Vg. The solid line is the fitted fm(yDc) parabolic curve from equation (4). Inset: the 
sideband reflection measurement scheme using a vr matching circuit. 

eling of the vr-matched LC resonator, the time varying Cg{t) leads to an observable frequency 
modulation of the tank circuit down to 6C ~ 0.1 aF. In our RF reflection measurement, the 
voltage applied on the input of the tt circuit is 



V(t) = Vdc + Vac cos(umt) + Vlc cos{u Lct). 



(1) 



Here Vdc is the constant voltage bias on the graphene capacitor. Vac actuates the membrane 
at ujm, and Vlc marks the probe RF signal at ojlc- The impedance modulation by the 
capacitance leads to side peaks V± in the reflected power at uJic^^m- In our case, Ci = 1.5 
pF, 6*2 ~ 50 fF, L ^ 15 nH, and the fitting of a measured tank resonance at fic = 5.89 
GHz yields Q ~ 80 corresponding an overall resistance oi R ^ 6 Q. Using circuit analysis 
with these parameters and the parallel plate approximation for the graphene capacitor, we 



obtain 

^ . (4.0 X lO'm- A)4^ = (4.0 X lO'm- A)|£ (2) 

for the relation between the graphene vibration amphtude and the sideband voltage. Here 
6D = -^ Jq A{y)dy is the average deflection representing the mechanical vibration ampli- 
tude and we have assumed that ^ is independent of position along the resonator. In our 
experiments, 6D <^ D^q, SC <^ Ceq- 

The electrical force acting on the resonator is 

m = ~\D=D.,v\t). (3) 

When \Vdc\ is small, Deq and the tension in the graphene do not change much with Vdc- 
Without tension the effective spring constant becomes KeffiVDc) = Kq — ^§^\D=D^qV^(., 
where Kq is the intrinsic spring constant of the graphene resonator. Consequently, there is 
a parabolic bias dependence given by 

UVDc)=fmOa-lVEc), (4) 

where 7 = j]^§^\D=Deg is a constant governed by the geometry of the resonator as well as 
the non-idealities of the graphene sheet, e.g. rippling and edges. 

The measurements were carried out in a high vacuum chamber dipped into liquid Helium 
(4.2 K). The actuator signal was taken from the internal reference of SR844 RF lock-in (up 
to 200 MHz). The RF carrier signal at ulc was injected through a circulator located on 
top of the cryostat to the sample, and the reflected sidebands were amplifled using Miteq 
low noise ampliflers (band 4-8 GHz). After down- mixing to the actuator frequency, the 
sideband amplitude was recorded by the RF lock-in using a time constant nock-in = 30 ms. 
Fig. |2] displays the measured sideband voltage as a function of the actuator frequency / 
and the DC bias voltage Vdc, with the highest mechanical resonance frequency around 178 
MHz obtained on sample #1 with length L = 0.7 fun and width W = 1 /zm. The concave 
parabolic fm dependence on Vjjc flts well Eq. (4). 

From Eqs. (1) and (3), the driving force at cOm is proportional to ^\j:,=Deq ' ^dc " Vac- 
In our experiments, nonlinear resonance behavior was found on all samples at large drives. 
Fig. [3Fa) shows a typical set of resonance curves displaying hardening Duffing behavior, 
when we change Vac at constant Vdc- Contrary to other samples, softening Duffing effect 
was observed on sample 7^1, as shown in the inset of Fig. [3ta). The hardening/softening 



behaviors correspond to the sign of the coefficient 03 of the restoring force aiu{t) + asu^{t) 
in the Duffing equation, where u{t) is the vibration displacement. Due to the competi- 
tion between elastic and capacitive mechanisms [M], a^ can be either positive or negative, 
resulting in the observed hardening and softening Duffing behaviors, respectively. The de- 
pendence of ai on the capacitance derivatives and intrinsic elastic parameters is different 
from that of as, and thus, there is no clear cut connection between fm{Vg) and harden- 
ing/softening Duffing behavior. However, the actual criterion for the cross-over between 
hardening and softening Duffing regimes is still unclear to us. From Fig. [Sta), with 
L = 1.5 /im, W = 2 fiin and D^g ^ 500 nm for sample #2, we can deduce the critical 
vibration amplitude 6Dhyst when the hysteresis emerges. For the red resonance curve in Fig. 
|3[a), we have V±/Vlc — 9.5 x 10^^ which according to Eq. (2), yields for the critical average 
vibration amplitude SDhyst ~ 20 pni. Using the approximations A{y) = ~'^°'^y^yl ' Amax 
and bD = j^ Jq A{y)dy [31] for the normalized basic mode shape, we obtain A^ax = 26D, 
where A^ax = max{A{y)) is the antinode amplitude of the mechanical resonator with length 
L. So the measured critical amplitude for the real graphene resonator is about 40 pm, which 
agrees quite well with the theoretical prediction of 50 pm [35]. On the same sample, the 
height of the minimum detectable resonance peak is about 4% of the hysteresis onset ampli- 
tude, which corresponds to 1.6 pm in the total vibration resolution. Taking into account the 
bandwidth of our measurement system B = l/riock-in ~ 33 Hz, we get for the sensitivity of 
our measurement a/S^ = 1.6 pm/v5 ^ 0.3 pm/vHz, which is ~ 10^ larger than the value 
of ~ fm/A/Hz achieved in the most sensitive detection schemes using SET or SSET [^I5B]. 
For the fundamental mode of sample #2, where fm ~ 57 MHz and Q ~ 1400, the estimated 
RMS amplitude of its thermal motion at 4K is A'^s = \/ m — ^fff )2 ~ ^-^ P'^- Using the 
same mode shape A{y), we have Af^, = [Uo '''"'T'^'^^ A±xdy]'^' = ^1^1.- There- 

Ath I 

fore, a sensitivity of ,"""" ^ 35 fm/vHz, about one order of magnitude smaller than our 
present a/S^, is needed to observe the thermal motion. Further improvements on sensitivity 
can be made, e.g., by using a cold amplifier or decreasing the gap between the graphene and 
the gate. 

Different driving forces usually result in different Q values for mechanical resonators 
[I2| [25] . In order to single out factors affecting the Q value, it is necessary to keep the 
driving force as fixed as possible. For this purpose, we carried out measurements where Vac 
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FIG. 3. (a) Hardening Duffing behavior sample #2 around 56.3 MHz (with Vdc = ~30 V and 
drive generator power stepped from -50 dBm to -30 dBm). Inset: softening Duffing effect addressed 
on sample #1 around 177 MHz (with V^ = 30 V and generator power increasing from -5 dBm to 15 
dBm) . (b) A series of resonance curves of a graphene mechanical resonator under constant driving 
and the corresponding QiVoc) dependence which fits Eq. (5), the displacement current dissipation 
model (red dashed line). 
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was changed according to Vdc, keeping the product Vdc ■ Vav constant; the prefactor in 
the drive ^\D=Deq changed only by 1% according to our estimates. Fig. k5[b) shows the 
resonance curves of sample #2 under a constant driving force in the linear regime. The Q 
values at different Vdc obtained by Lorentzian fits are also plotted in Fig. p[b). The data 
shows a decreasing Q{Vdc) with increasing \Vdc'\ while the amplitude h of the response 
remain constant. As Q ■ h increases, this means that either the mode shape changes, or 
there is a slight increase in ^\D=Deg- 

Typically, tension enhances the mechanical quality factor. In graphene stretched by 
electrical forces, displacement currents may induce considerable dissipation due to the large 
resistance of graphene. The displacement current generated by Vdc and time varying capac- 
itance is given by IdiVDc) = ^mVDc^C ■ siniumt) . Since the resistance of gated graphene be- 
haves as RciVDc) = [5'o(l + /3| Vdc|)]~\ the energy dissipated per period by the displacement 
current becomes i?e(^Dc) ~ /o ™ -^J-^C^^ = ^^^^^^^^"^ — " i+my — i' ^^'^^ 9^ ^^ ^^^ minimum 
conductance and (3 parametrizes the increase in conductance with Vdc- In Fig- [3|^b), Um 
does not change much with Vdci ciiid the total energy stored in the resonator Et can be taken 
as a constant due to the equal heights of the resonance peaks. So with a Voc-independent 
energy loss Em per period, the overall Q value can be written as Q{Vdc) = ^ +e^{v — T' '-'^ 

Q{Vdc) Qm 1 + /3|Vdc|' ^^ 

where Q^ = Q{0) = 27iEt/Em and a = ^^^^ = ^^^'.JZ2' ^^^ constants. With go = 
23 /iS and /3 = 0.01, Eq. (5) fits well with the experimental data, as shown by the red 
dashed curve in Fig. p[b). For this sample, we obtain wig// = 2.2 x 10~^^kg by fitting 
fmiVDc) with Eq. (4). For the effective mass density we obtain pe// = Q.Tpo, where 
Po = 7.6 X lO^^^kg/pm^ is the ideal mass density of monolayer graphene. Note that the 
suspended graphene conductance is rather insensitive to the gate voltage (i.e., (3 is small) 
due to the absence of a dielectric layer between the gate and graphene. As a consequence, 
there is no need in Eq. (5) to take into account the Dirac point shift, which would bring 
only an insignificant asymmetry to the QiVDc) curve. The large pejf is probably a result 
of unpredictable adsorbates, while the small go may arise due to bad contact resistance in 
addition to adsorbates. Improvement on these numbers can be expected by annealing |23] . 
However, our method is easily generalized, for graphene mechanical resonators suspended 
on rigid target supports by using a similar non-flip assembly. 



Complementary to the down-shifting behavior of Eq. (4), the deformation-induced ten- 
sion in graphene leads to up-shifting tendency at higher |Vdc|- Iii Fig- R we present data 
from a few-layer graphene sample #3 with L = W = 1 /im, where the "W"-shape fmiync) 
curve shows a clear transition from the electric-force dominated regime to the tension dom- 
inated behavior. In the high \Vdc\ regime, fmiVoc) fits the fm oc V^^ model [HI 137] , 
indicating a negligible intrinsic tension in the transferred graphene of this sample. In the 
intermediate |Vdc| regime, fm is less dependent on Vnc and the resonance response curves 
become wider and asymmetric, implying some extra dissipation and nonlinear mechanics. 
Analogous experimental signatures in this transition regime have been found in a recent 
work on nanowire resonators [3H], where the authors attribute the widening of resonance to 
the mixing of different modes. In our case of graphene, where no other modes were involved 
in the adjacent frequencies, a convincing theoretical model remains to be worked out. In 
the low \Vdc\ regime, fm{VDc) fits well with the parabolic model of Eq. (4), which gives 
an effective mass wig// ~ 5 x 10^^^ kg for this sample, corresponding to pe// ~ 0.65po when 
Q]^\D=Deq is calculated from the parallel plate capacitance model. Although it is known 
that the effective mass of a clamped-clamped mechanical resonator can be smaller than the 
real mass of the suspended part [21], we attribute the small determined mass to edge effects 



that increase the actual value of |^ than the calculated coefficient. In the case of clamped- 
clamped graphene, the unpredictable adsorbates and free edges make it more difficult to 
get an ideal effective mass due to the extremely small mass and thickness of graphene it- 
self. Hence, we conclude that, adsorbates and edges must be carefully considered in both 
sample-making and modeling to determine the resonator mass correctly. 

In conclusion, we have developed a micro-stamp technique that can be used to transfer 
and assemble suspended graphene into nanoscale mechanical resonators. This new technique 
enabled us to take advantage of a localized gate in close proximity to the sample, which 
improved the sensitivity of the RF-cavity readout scheme by enhancing the capacitance 
ratio between the graphene device and the parasitic contributions. We found Duffing effects 
of different sign on graphene samples, and measured the critical vibration amplitude at 
which the Duffing hysteresis emerges. We observed "W"-shaped, nonmonotonic fmiVnc) 
curve on a few-layer graphene mechanical resonator, demonstrating a continuous transition 
from electrically dominated low \Vdc'\ regime to tension-dominated high \Vdc\ regime. In 
our experiments on a series of graphene samples, we obtained resonance frequencies up to 
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FIG. 4. The mechanical resonance frequency fm as a function of DC gate voltage Vdc for sample 
7^3. The amplitude of the sideband V± is given by the color scale in Volts. The non-monotonic, 
"W" -shaped curve shows a transition from the capacitor-dominated (low |Vdc|) regime to the 
tension-dominated (high |Vdc|) regime. The other two straight traces with little Vdc dependence 
are from the gold clamps with hanging ends like cantilevers. Inset: the two red dashed lines are 
the two different fitting curves for low/high |Vdc| regimes, respectively. 



178 MHz, a sensitivity of about 0.3 pm/vHz, and an effective resonator mass down to 
~ 10^^^ kg. The achieved combination of high frequency, high sensitivity, and low mass 
show that our micro-stamp technique and RF reflection measurement scheme hold promise 
for building novel, sensitive NEMS structures. By using symmetrized electrical cavities |39]. 
the sensitivity of our scheme can be further increased and the quantum limit is within reach 
with graphene mechanical resonators. 
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